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ABSTRACT 

We present follow-up observation^] and analysis of the recently discovered short period low-mass 
eclipsing binary, SDSS J001641-000925. With an orbital period of 0.19856 days, this system has one 
of the shortest known periods for an M dwarf binary system. Medium-resolution spectroscopy and 
multi-band photometry for the system are presented. Markov chain Monte Carlo modeling of the 
light curves and radial velocities yields estimated masses for the stars of Mi = 0.54 ± 0.07 M Q and 
M 2 = 0.34 ± 0.04 M Q , and radii of Ri = 0.68 ± .03 R© and R 2 = 0.58 ± 0.03 R respectively. This 
solution places both components above the critical Roche overfill limit, providing strong evidence 
that SDSS J001641-000925 is the first verified M-dwarf contact binary system. Within the follow-up 
spectroscopy we find signatures of non-solid body rotation velocities, which we interpret as evidence for 
mass transfer or loss within the system. In addition, our photometry samples the system over 9 years, 
and we find strong evidence for period decay at the rate of P ~8 s yr _1 . Both of these signatures raise 
the intriguing possibility that the system is in over-contact, and actively losing angular momentum, 
likely through mass loss. This places SDSS J001641-000925 as not just the first M-dwarf over-contact 
binary, but one of the few systems of any spectral type known to be actively undergoing coalescence. 
Further study SDSS J001641-000925 is on-going to verify the nature of the system, which may prove 
to be a unique astrophysical laboratory. 

Subject headings: binaries: eclipsing - binaries: spectroscopic - stars: low-mass 



1. INTRODUCTION 

M dwarfs, with masses between ^0.08 and 0.6 M , 
are the most numerous and faintest main sequence stars 
in the Galaxy. The high number density and very long 
main sequence lifetimes of M dwarfs allow them to be 
used as det ailed tracers of the nearby g alactic stellar 
population (jBochanski et al.|[2007al 120101) . These stars 
are also well known for their high magnetic activity lev- 
els, characterized by frequent violent flaring out bursts 
(e.g. iKowalski etaD 120091; iDayenport et all 120121), an d 
Ha emission (jWest et alll2008l 120111: Iffilton et alfeoiol ). 

Eclipsing binary systems provide us with the valuable 
ability to directly measure the mass es and radii of t he tw o 
component stars. As discussed in iBecker et al.l (| 2 1 If ) . 
the sample of known binary systems composed of two M 
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dwarfs is very small. Their intrinsically low luminosity 
and small radii conspire to make M dwarf binaries dif- 
ficult to detect in variability surveys. The fraction of 
M dwarfs in binary systems is a lso much lower than for 
higher mass stars. ILadal ()2006D has estimated that as 
few as ~10% of M dwarfs in the Galactic field may be in 
binary (or higher number) configurations. 

The M dwarf binary and single-star systems that are 
well studied indicate that stellar structur e models do 
not accurately pre dict the observed radii (jRibasI 120061 : 
iMorales et alJl2010h . Radii are often observed to be ~10- 
20% larger than predicted for a given mass. This may be 
due to limitations of the stellar models used, or to a fun- 
damental difference in the interior structure of M dwarfs. 
Magnetic activity is frequently identified as a potential 
culprit, causing the radii to be enlarged due to the inhi- 
bition of co nvection or decreased heat flux from magneti c 
spots (e.g. iChabrier et all [20071 iBochanski et all 1201 ID . 
Long period M dwarf binaries (P > 10 days) are ideal 
systems for radii determinations, as they are no t influ- 
enced by spin-orbit coupling or tidal distortion (|Mazehl 
2008). The above discrepancies will only be resolved by 
the addition of many accurately measured radii to the 
known sample over a wide range in mass and period. 

There also appears to be an absolute lower-limit on the 
orbital period for contact binary star systems, though 
no consensus about its origin has been reached. This 
short period limit at ^0.22 days may be due to both 
components being fully convective (jRucins O ll992l ). At 
shorter periods, these stars cannot adhere to their mass- 
radius relationship while simultaneously stayi ng within 
the Ro che equipotential geometry. However, as iRucinskH 
(1992) points out, this does not appear to be the primary 
effect at work in creating the sharp 0.22 day period limit. 
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Figure 1. Phased light curves in the five SDSS filter bandpasscs 
for SDSS J001641-000925. Photometric errors are shown for all 
data. Blue points are from the Stripe 82 catalog, black arc follow- 
up data from the NMSU 1-m in griz bands. Red data in the i-band 
were obtained with the Agile camera at the APO 3.5-m. 

Alternatively, iJiang et al.l (|2012f ) have shown that unsta- 
ble mass transfer for low mass contact binaries may result 
in rapid coalescence. The observed period limit would 
thus be a result of the very short lifetimes for such sys- 
tems. Models of angular momentum loss in binary sys- 
tems, however, indicate that even given favorable mass 
ratios and initial orbital periods, 6-13 Gyr timescales are 
typically required for l ow mass binari es to obtain periods 
as short as 0.22 days (jStepienl 1200 6b). The observed pe- 
riod limit would then be a consequence of the finite age 
of the binary population. 

Until recently, the shortest known M dwarf bi- 
nary system period in the literature was BW3 
V38 flMaceroni fc Rucinskilll997tlMaceroni fe Montalbanl 
I2004D . a detached binary composed of two main se- 
quence M3 dwarfs w i th an orbital period of 0.1984 days. 
iRucinski fc Pribullal (|2008| ) found a similarly short pe- 
riod system, GSC 01387-00475, with a period of 0.2178 
days, but composed of higher mass K3 and K5 stars and 
i n a contact configuration. 

iDimitrov fc K iurkchicva (2010) reported the shortest 
period M dwarf binary yet characterized, GSC 2314-0530 
with a 0.192636 day period, and stellar masses of 0.51 
Mq and 0.26 M . This system has an orbital period 
only ^8 minutes shorter th a n tha t of our target, SDSS 
J001641-000925. iNefs et al.l ([20121 ) have recently discov- 
ered four binary systems with M dwarf colors that have 
periods significantly below the ~ 0.22 day period limit, 
providing a major challenge to our understanding of the 
formation and lifetimes of such systems. The rarity of 
short period M dwarf binaries, and the lack of known 
contact M dwarf binary systems, indicates that they are 
difficult to produce, are so unstable that they are rapidly 
destroyed, or are still very difficult to detect. 

A recent stud y of photometric variability by 
IBecker et al.l (|2011l he reafter Bll) usin g the Sloan Digi- 
tal Sky Survey (SPSS. l^rFeTall2000Tl revealed several 
new candidate M dwarf - M dwarf eclipsing binary sys- 
tems. Bll analyzed ^4.8 million light curves and found 
207 periodically variable objects, and a sharp drop in the 
number of systems with periods less than 0.22 days. Of 



these, two candidate eclipsing binary systems with sinu- 
soidal light curves and M dwarf ugriz colors were found 
to have very short periods: SDSS J20001 1.19+003806.5 
at 0.14552 days, and SDSS J001641.03-000925.2 at 
0.198561 days. Both of these objects warranted spec- 
troscopic and photometric follow up to determine their 
fundamental parameters, and to confirm their status as 
short period eclipsing binaries. 

In this paper we describe our detailed photometric 
and spectroscopic investigation of the brighter of these 
two objects, the 0.198561 day periodic variable SDSS 
J001641-000925. In §2 we describe the photometric and 
spectroscopic data used in this study. We model basic 
parameters in §3, and analyze the peculiarities of the 
binary system in §4. We discuss the context of this over- 
contact binary in §5. 

2. OBSERVATIONS 

SDSS J001641-000925 was discovered as a peri- 
odic variable object by Bll using time-domain data 
from the SPSS "Stripe 82" region (jlvezic et al.l 120071: 
iSesar et al.l 120071 ) . This unresolved binary system has 
an SDSS ObjID=587731 186205720773, and SpecOb- 
jID=314914764310446080. A period of 0.198561 days 
was determined from 61 epochs of griz photometry, span- 
ning more than 8 years of observations, using the Su - 
perSmoother algorithm (Friedman 1984; Rcimann 1994). 
This algorithm does not assume a light curve shape, 
and instead works to find the period that minimizes 
the phase folded light curve dispersion. SDSS J001641- 
000925 has a median r-band magnitude of 16.39. A pho- 
tometric spectral type of Ml was found using the me- 
dian riz-bands and th e covariance matrix technique from 
iKowalski et al.l (|2009l ). while a type of M0 was assigned 
using the "Hamm er" spectral-typing analysis package 
(jCovev et al.l 120071) . which compares spectral data to a 
set of standard templates. The SDSS spectrum did not 
show any clear signatures of binarity, such as spectral 
line separ ation or large deviation from the M0 spectral 
template (Bochans ki et al.ll2007bO . In the following sec- 
tion we describe follow-up data obtained to characterize 
this short period object. 

2.1. Photometry 

Figure[T]contains all of the photometry for this system. 
In addition to the 61 epochs of SDSS Stripe 82 photom- 
etry, we obtained follow-up i-band photometry from the 
Agile high speed imaging camera on the Astrophysical 
Research Consortium 3.5-m telescope at Apache Point 
Observatory (APO), and griz-hand photometry from the 
New Mexico State University (NMSU) 1-m. 

The SDSS photometry comes from the "Stripe 82" re- 
gion, a narrow strip located along the celestial equator 
(20 h 24 m < aj 2 ooo < 04 h 08 m , |<5 J20 oo| < 1-266°) that has 
been repeatedly observed throughout the course of the 
SDSS. This catalog produced light curves for almost 8 
million point sources with between ^20 and ~80 epochs. 
These data have been used for a va riety of investigations , 
including stud ies of supernoyae dFrieman et al l [2008) , 
stellar flares (Kowalsk i et al.l I2009D . and characteriz- 
ing m ultiwavelength point source variability (|Sesar et al.l 
12007th SDSS J001641-000925 was observed 61 times in 
all five photometric bands. 
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Figure 2. One of the nine MagE spectra for SDSS J001641.03- 
000925.2, taken at a phase of = 0.236. The spectrum has been 
smoothed using a boxcar kernel of 7 pixels. 

The Agile camera observed the system for ^45% of the 
orbital period (2.2 hours) on 2008 October 30, with a ca- 
dence of 90 seconds, and a median photometric error of 
0.004 magnitudes. The ^0.2 magnitude sinusoidal vari- 
ation in the Agile light curve appeared in phase with the 
sparsely sampled Stripe 82 data (see Figure Q}. 

Multi-band follow- up was provided by th e automated 
NMSU 1-m telescope (jHoltzman et alJl2010f >. These griz 
data were taken over a series of 16 half-nights in 2010 
from November 13 to December 2. The observing strat- 
egy and exposure times varied between the 16 nights. 
During two nights we imaged the binary with consec- 
utive griz filter exposures. Single filters were used for 
each of the remaining 14 nights, broken down as follows: 
two nights of g, seven nights of r, three nights of i, and 
two nights of z-band. These phased light curves matched 
the well defined sinusoidal trends seen in the griz-bands 
from the SDSS data, with no major deviations e.g. from 
flares or star spots. 

2.2. Spectroscopy 

SDSS spectra are the result of co-adding three to 
six in dividual frames, e ach with 9-15 minute exposure 
times (|York et all [20001. The time-domain SDSS spec- 
tra have previously been examined t o search for emis- 
sion line variabilit y in M dwarfs (e .g. iHilton et al.ll2010t 
iKruse et al.ll201Ct iBell et al.ll2012ft . and quasar variabil- 
itv (jVanden Berk et al.l 120041 ). We obtained three time 
resolved spectra from 2002 October 12 using the SDSS 
Data Archive ServerQ No spectral line splitting was ap- 
parent from the R ~ 1800 resolution SDSS spectrum for 
this system. A single Ha emission line was seen in each 
frame, but no radial velocity shift was resolved for any 
spectral feature. 

Initial follow up spectroscopy was provided by the 6.5- 
m Magellan Clay telescope at Las Campanas Observa- 
tory, using the R ~ 4100 resolution Magellan Echellette 
Spectrograph (MagE) on 2008 October 25. We obtained 
nine spectra with exposure times ranging from 600 to 
900 seconds. These data were reduced using the MagE 
Spectral Extractor software package for IDL (MASE; 

7 http : //das . sdss . org 
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Figure 3. Nine MagE spectra showing clear splitting in the Ca I 
lines at 6102A and 6122A. The MagE spectra were taken in order 
from top to bottom, and cover 44% of the orbital phase, from 
</> = -0.15 to = 0.29. 
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Figure 4. Radial velocities as a function of orbital phase for both 
components of this system for the MagE and DIS data. The veloc- 
ities for the primary star are open black circles, and the secondary 
star are solid grey circles. The best-fit model, as generated by the 
PHOEBE software, is overlaid. 

iBochanski et al.ll2009D . and an example is shown in Fig- 
ure [2] The observations spanned ^44% of the orbital 
phase. 

Additional spectroscopy was obtained with the APO 
3.5-m Dual Imaging Spectrograph (DIS), using the 
R1200 grating with a resolution of R ~ 3500, on 2011 
October 20. These eight spectra had integration times 
of 900 seconds, and were timed to cover ~50% of the 
orbital phase that was not observed with MagE. 

3. MODELING THE SYSTEM 

Using the M dwar f spectral templates from 
Bochansk i et al.l (|2007bT ) we estimated the spectral 
type for both stars in this system. We scaled these 
template spectra, from types M0 to M3, by their approx - 
imate bolometric luminosities (R eid fc Hawleyl [2005), 
and created a normalized binary template spectrum for 
each of the 16 possible combinations of stars (M0+M0, 
M0+M1, ... M3+M3). The average best fit template 
combination for the DIS and MagE data was M0+M1, 
which we adopted to seed the initial binary model. 

Velocities were measured for both stellar components 
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Figure 5. Top: Probability distributions of the radii normalized 
by the semi-major axis, as a function of the mass ratio, for the 
primary star (blue contours) and secondary star (red contours) 
explored by the MCMC sampler using the Unconstrained mode in 
PHOEBE. The solid blue and r ed lines show the effective Roche 
lobe radii from lEggletonl ((1983) for the primary and secondary 
stars, respectively. The effective radii for both stars are larger 
than the critical Roche lobe radii, indicating an over-contact binary 
scenario. Bottom: A three dimensional model for SDSS J001641- 
000925 at a phase of <j> = 0.25 from the PHOEBE software using 
our best-fit values, with the primary star (left), secondary star 
(right), center of mass (cross), and the critical Roche equipotential 
(solid line). 

by cross correl ating the active MO an d Ml spectral 
templates from .Bochanski e t al.l (|2007bf) to each of the 
MagE and DIS spectra using a two-dimensional cross 
correlation code empl oying the TODCOR algorithm 
(|Zucker fc Mazehlll994P i. Figure [3] presents three absorp- 
tion lines in the MagE data that show clear radial ve- 
locity splitting, indicative of a binary system. Veloci- 
ties were also extracted from the spectra usin g the IDL 
broadening function code from Rucinski (2002), and were 
consistent with the cross correlation results within the 
~15 km s _1 TODCOR uncertainties. The radial veloc- 
ity curve for the MagE and DIS data is shown in Figure 
H] as a function of orbital phase, with the best-fit binary 
model (described below) overlaid. 

3.1. Initial Binary Model 

In order to combine the many sources of data, we used 
the latest st able build of the PH OEBFfl modeling soft- 
ware ( v0.32;|P rsa fc Zwi tter 2005). This implementation 
of the I Wilson fc Devinnevl (|197lf) algorithm is an indus- 
try standard in modeling binaries, and is actively being 

8 http://phoebe.fiz.uni-lj.si/ 



Figure 6. Same as Figure [5] but for the W UMa Overcontact 
mode in PHOEBE. 

deve loped and improv ed for use in modern binary studies 
(e.g. iPrsa et aLll2011[) from missions such as the Kepler 
space telescope (jBorucki et al.|[2010l ). 

Due to the very short orbital period of the system we 
initially used the ' 'Undetermined' ' binary configura- 
tion in the PHOEBE software, which allowed solutions 
wherein both stars could potentially overflow their Roche 
lobes. We estimated the initial system parameters as- 
suming a M0+M1 binary pair, as found by our spectro- 
scopic fit described above, with masses of 0.6 and 0.5M Q , 
temperatures of 3800 and 3600 K, and an orbital sepa- 
ration of 1.4 Rq. 

The binary parameters were converged manually using 
the "Differential Corrections" solution mode. We started 
by fitting only the photometry, assuming 70 degrees in- 
clination and a circular orbit. The light curves for the 
griz bands were fit simultaneously, and the effective tem- 
peratures and luminosities of the stars were estimated. 
The radial velocities from the MagE and DIS data were 
then included in the model, and the solution recomputed. 
The orbital parameters were then allowed to vary, and 
the solution was recomputed again using both the ve- 
locities and photometry as constraints. Eccentricity re- 
mained fixed at 0. The orbital period from Bll for the 
system was adopted, with no time derivative or phase 
shift in the period. A logarithmic limb darkening pro- 
file, included with PHOEBE, was used. No surface spots 
were included in our binary solution. 

3.2. MCMC 

The result of our initial binary solution indicated that 
the system was in a contact configuration, with the pri- 
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Figure 7. Probability distributions for the masses, radii, and the semi-major axis for SDSS J001641-000925 as explored by our MCMC 
sampler using the W UMa Overcontact mode in PHOEBE. 

mary star overflowing its Roche lobe. This would make 
SDSS J00f641-000925 the first contact M dwarf binary 
system yet characterized. To determine the reliability 
of the contact solution, we employed a Markov Chain 
Monte Carlo (MCMC) parameter space search method. 

We used the Python-base d MCMC sampler emcee 
(jForeman-Mackev et alJl2012T ) to adjust the free param- 
eters Ti, T2, a, i , and t he potentials for the two 
stars Qi and f2 2 {Kopal 1959). Initial uncertainties for 
our MCMC runs were determined using the error pre- 
scription included in the PHOEBE v0.32 user manual 
for our initial binary model parameters. We then used 
the command line PHOEBE "Scripter" to compute the 



X", masses, and radii for each iteration (step). 

A total of 16 parameter space "Walkers" were used, 
and multi-threading with seven computer cores was im- 
plemented using the Python multiprocessing module. 
The first 1,000 steps were discarded, to allow the sam- 
pler an initial exploration of the parameter space. The 
MCMC was run for an additional 20,000 steps for each 
Walker, for a total of 320,000 iterations. 

We ran our MCMC search using both the 
"Overcontact Binary of the W UMa type" mode 
and "Unconstrained" mode in PHOEBE, as the con- 
straints for each mode can affect the resu lting parameter 
estimates, especially of the radii (e.g. iBonanosl l2009t 
iCaballero-Garcia et al.l |201Q[ ) . The best-fit results for 
the free parameters in both modes were the same 
within the uncertainties. The secondary star radius, 



however, was very different between the modes. In 
Figure [5] we show the radii versus mass ratio probability 
distributions, and a three dimensional illustration of the 
system at a phase of <fi — 0.25, with the Unconstrained 
mode MCMC runs. Figure U shows the MCMC results 
using the W UMa Ove rcontact mode. We note t hat the 
Over contact model (Wils on fe Devinnevl Il97lt iWilsonl 
|2001[ ) can converge on a phys i cally incorrect solution 
(e.g. see iPribulla fc Rucinskil 120081 ). particularly for 
systems with mass transfer or mass loss as discussed in 
§4.2. Both mo des, however, indicated that t he system 
is in contact (jRucinskil 119971; IWilsonl 120011) . and we 
present results for the Overcontact W UMa mode as 
our preferred final model. The stars were found to have 
masses of Mi = 0.54 ± 0.07 M and M 2 = 0.34 ± 0.04 
M(T), which correspon d to spectral types Ml and M3 
(|Reid fc Hawlevll2005h . and radii of Ri = 0.68 ± .03 
R and R 2 = 0.58 ± 0.03 R . A grid of the resulting 
probability distributions for the masses, radii, and 
semi-major axis is shown in Figure [7] 

We computed the correlation lengths for each of our 
free parameters using the autocorrelation function. The 
effective chain length is defined as the actual length of the 
MCMC chain after the 1,000 step "burn-in", divided by 
the correla tion length, and mus t be >> 1 for meaning- 
ful results (jTegmark et al. 2004). The shortest effective 
chain length in our study was for the inclination param- 
eter, with a value of 7,106. All other free parameters 
had effective chain lengths over 9,000, indicating that 
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Figure 8. The differential SDSS g-band, (g — r), and (r — i) colors 
as a function of orbital phase, along with the PHOEBE model fit 
(red). 



Table 1 

The system parameters and uncertainties 
for SDSS J001641-000925 determined by 
our MCMC runs. 



Quantity 


Value 


Uncertainty 


P (day) 


0.198561 


lxl0~ 6 


t (day) 


0.002405 


lxl0~ 6 


Mi (M ) 


0.54 


0.07 


M 2 (M ) 


0.34 


0.04 


Ri R Q ) 


0.68 


0.03 


R 2 (R ) 


0.58 


0.03 


T cff ,i (K) 


4342. 


475. 


T eff 9 (K) 


3889. 


579. 


a (R ) 


1.37 


0.06 


q (Ma/Mi) 


0.62 


0.01 


i (deg) 


53.3 


1.1 


v aya (km s _1 ) 


-5.6 


3.9 



our parameter space sampling reliably converged on the 
correct solution, and that the resulting errors were real- 
istic. We also tested the convergence of our MCMC runs 
with the Gelman-Rubin R statistic, as detailed in §3.3 
of iFordl ((2006). All parameters converged to R < 1.001, 
indicating our MCMC effectively sampled the parameter 
space. 

The final radial velocity solution is shown in Figure 
HI and the SDSS phase curves computed from the final 
PHOEBE model are shown in Figure [8] The median 
model values for the system parameters are given in Ta- 
ble [TJ Errors for each parameter were determined using 
the standard deviation of the probability distributions. 

4. RESULTS 

4.1. Binary System Parameters 

The illustration of our binary system in Figure[S]clearly 
indicates that both stars are filling their Roche lobes. We 
also show the probability distributions for the normal- 
ized radii of both stars. Overlaid is the critical Roche 
radius for both t he primary and secondary stars, cal- 
culated using the lEggletonl (|1983l ) approximation. The 
MCMC probability distribution for the W UMa model 
mode showed both stars to be larger than their critical 
Roche radius with 100% confidence. The primary star 
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Figure 9. The stellar radius as a func tion of mass fo r the c om- 
pilation of single and b inary stars fro m Lopez-Morales (20Q3), as 
well as the system from Becker ct al. (2008), in open black circles. 
Our PHOEBE solution for SDSS J00164 1-000925 is shown as filled 
blue squares. The solid black line is the Baraffc ct al. (1998) stel- 
lar evolution model, an d the dashed line is the empirical fit from 
IBavless fc Oroszl pOOlft . 

was 18% larger than the critical Roche radius, and the 
secondary star 25% larger. This raises the i ntriguing 
possi bility that the system is in over-contact (jRucinskil 
|l997j), in which case we would expect to see non-solid 
body velocities from outflowing material leaving the sys- 
tem via the L2 point. This would be a very efficient 
method of angular momentum loss for the binary, and 
could lead to rapid decay of its orbit. 

The masses and radii of SDSS J001641-000925 with 
respect to o ther low-mass ec l ipsing binaries are plotted 
in Figure [S IBavless fc Oros z (2006) provided an empiri- 
cal fit to the Mass-Radius relation using many low-mass 
binary systems with larger tha n predicted radii . Our 
system is much larger than the IBavless fc Oroszl (|2006[ ) 
prediction at this mass range, due to the tidal distortion 
effects of the over-contact solution. 

4.2. Ha Line Emission 

Ha emission from the system was present in every ob- 
served spectrum. This is not an uncommon property 
of M dwarfs in general, whose magnetic activity is fre- 
quently chara cterized by strong and variable Ha line 
emission (e.g. iWest et all 120111 : Bell et all 120121 ). How- 
ever, the fraction of magnetically active fiel d stars with 
spect ral types M2 or earlier is less than 10% (|West et al.l 
1201 If ). The strength and frequency of this emission de- 
creases as M dwarfs age, and the lifetime of magnetic 
activity f or early typ e M d warfs in the field is less than 
-2 Gyr (IWest et all 120081) . Ha emission is also seen 
from M d warfs in close binary systems w ith other dwarf 
stars (e. g.lDimitrov fc Kiurkchievall2010f) and with white 
dwarfs (jMaxted et al.l 119981: ISilvestri et all I2006D . The 
higher fraction of active early type M dwarfs in close 
binaries as compared to the field is likely due to these 
stars rotating more rapidly than their isolated counter- 
par ts, inducing stronger and longer lived mag netic activ- 
ity (|Morgan et alj|2012t ISilvestri et aHl2006f ). 

We measured the Ha equivalent width (EW) for each 
DIS, MagE, and SDSS spectra, shown in Figure flOk as 
a function of orbital phase. The strength of emission is 
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Figure 10. a) The equivalent width of the Ho line in each spec- 
trum as a function of orbital phase for the MagE (filled black cir- 
cles), DIS (open blue squares), and SDSS time-spectra (open red 
diamonds), b) The radial velocity of the Ha line for each spectrum, 
with the orbital radial velocity solution from Figure [4] overlaid for 
comparison, using the same color coding and symbols as above. 

nearly constant at ~lA, though a weak correlation in the 
line strength was found with orbital phase, with a higher 
Ha EW at conjuncture (<j> — and <fr — 0.5). 

The Ha line was observed to be broader than other 
features (e.g. the Ca I absorption lines in Figure [3]), and 
clear splitting was only seen in a few spectra. The feature 
did not contribute significantly to the cross correlation 
in §3, which was dominated by absorption features such 
as those shown in Figure [3] We measured the radial ve- 
locity of the Ha line, shown in Figure [TUb . by fitting two 
Gaussian profiles to the emission feature. The Ha line 
broadly followed the velocity trend of the more massive 
stellar component, possibly due to the primary domi- 
nating the total Ha flux. Line splitting was seen near 
quadrature (<j> — 0.25) in the higher resolution MagE 
spectra, and the Ha velocity profiles are consistent with 
emission from both stellar components in this system. 

No other evidence for accretion, flaring, or outflow in 
the form of nebular emission line features was observed. 
The line splitting seen in the MagE data at <j> = 0.25 indi- 
cates that both stars are magnetically active. We would 
conventionally assume the Ha emission comes from many 
distributed active regions on the stellar surfaces. How- 
ever, the broadened Ha emission, lower amplitude radial 
velocities, and intermittent line splitting suggests non- 
solid body velocities, and may indicate mass transfer be- 
tween the stars and mass loss, likely through the L2 point 
as a result of the over-contact configuration. The slight 
increases in Ha EW at conjuncture are reminiscent of 
those from the massive over-contact binary system RY 
Scuti, wh ich are attributed to exten ded Ha emission from 
outflows (|Grundstrom et all 120071 ) that appear to have 
brighter EW during eclipse (0 = and <fi = 0.5) when 
the system flux is lower. Higher resolution spectroscopic 
monitoring throughout the orbital period is necessary to 
accurately determine the geometry of the Ha emission 
regions, and the possible active region timescales. 
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Figure 11. The difference in orbital period between the Bll value 
and our measurements, using SDSS (blue) and NMSU (black) pho- 
tometry, as a function of time. Periods shorter than the Bll value 
yield negative AP. The first SDSS epoch, at At = 0, had a MJD 
= 51075.316. Horizontal bars give the range of data for each SDSS 
subset, with each point placed at the mean epoch for the subset. 
Error bars were determined using a bootstrap method, described 
in the text. 



4.3. Period Decay 

The observed period of a close binary system should 
decrease over time, as the orbit decays and the system 
loses angular momentum from effects such as magnetic 
braking, tidal friction, mass transfer, mass loss, and in 
the most compact cases gravitational waves. Changes in 
the orbital period have been obser ved for several hun- 
dred W UMa-type binary systems. iPilecki et all (2007) 
measured period decreases for ~30 short period sys- 
tems, finding typical amplitudes of P = =£■ ~ — 1 
sec year - 1 iKubiak et all d2006ft searched the OGLE 
database (|Udalski et all 119921) for period changes from 
more than 500 short period binaries, finding a maximum 
amplitude of P = — 0.4 sec year -1 . 

The SDSS Stripe 82 data were very sparsely sampled in 
time, and individual times of light curve minima were not 
observed. As such, we were unable to make a traditional 
Observed - Computed (O-C) diagram. Instead we broke 
the SDSS data into five subsets, spanning different ranges 
in time, each having at least 20 epochs. The NMSU 1-m 
photometry were obtained in a single observing season, 
with many full orbits observed, and were used as a sixth 
subset. 

To search for changes in the orbital period of SDSS 
J001641-000925, we measured the period for each time 
interval in our data with a phase dispersion minimiza- 
tion c ode, using the Su perSmoother alg orithm ()Friedmanl 
119841 ) as described in iReimaiml (|1994|) . We folded the 
each subset at periods ranging ±5 minutes around the 
Bll value, in steps of 0.1 seconds. 

In Figure fTTI we show the difference in period from the 
Bll value for each subset of r-band Stripe 82 and NMSU 
data. The time span for each Stripe 82 subset are shown 
as horizontal bars. A consistently decreasing period was 
found, with a total amplitude of ^75 seconds, spanning 
the 9 years of observations, yielding P ~ —8 sec year -1 . 
A similar app roach was used to me asure period decay for 
V1309 Sco in iTvlenda et all poTl . 



8 



Davenport et al. 



We recomputed the period determination for each data 
subset using a bootstrap method, leaving one epoch out 
each time. Errors shown in Figure [TT] are the standard 
deviation of the bootstrap runs. The Stripe 82 period re- 
sults are somewhat dependent on the time range chosen, 
as shown in Figure 1111 The independent period mea- 
surement from the NMSU photometry, however, provides 
strong evidence of orbital decay for this system, though 
we note the caveat that this is dominated by the NMSU 
photometry. 

5. DISCUSSION 

We have presented follow up spectroscopic and photo- 
metric observations of a low-mass binary, SDSS J001641- 
000925, whose orb ital period is below the theoretical 
short period limit ([Rucinskil 119921 ) . The components of 
this short period system have masses typical of spectral 
types of Ml and M3, and are in an over-contact configu- 
ration. Both stars have radii much larger than predicted 
by stellar isochrone models. These stars are the most 
"oversized" of any yet determined at this mass due to 
the over-contact configuration. 

The strength of the Ha line at ~lA was not partic- 
ularly high with respect to active stars of this spec- 
tral type. Higher Ha flux is typically accompanied by 
strong magnetic activity and flaring for early type M 
dwarfs; these are features that diminish as the st ars age 
and lose angular momentum ([West et al.l 12008ft . The 
SDSS spectroscopic M dwarf sample, c ontaining ~70,000 
spectroscopically confirmed M dwarfs ([West et al.ll201lft 
shows approximately half of stars with spectral types 
M0-M3 that are classified as "active" have EW(Ha) > 
1.25A. Thus the Ha EW for this system is not usu- 
ally high, nor does it appear particularly active. In 
fact, no signs of flaring were detected in more than 
20 nights of observing. This is in strong contrast 
to the si milar mass semi-detached bina ry system ob- 
served bv lDimitrov fc Kiurkchieval {20 lO), which showed 
six flares in a c o mpar able amount of observing time. 
iDavenport et al.l ([20 12) found that field stars in this 
spectral type range displayed photometric flares with am- 
plitudes of Am > 1 mag once every ~2 days on average. 

However, the slight increase in Ha EW during eclipse, 
as well as the broadened Ha line profile with possibly 
non-solid body rotation velocities, together suggest that 
the emission feature may be the result of material be- 
ing expelled from the outer Lagrange points. This would 
result in a complex Ha line, with phase-dependent emis- 
sion lines contributed from both stars, as well as possible 
P Cygni-like profiles from the outflowing material. Addi- 
tional medium to high resolution spectroscopic monitor- 
ing will allow us to differentiate Ha emission due to active 
surface regions from that of any ne bular emission sur- 
round ing the system, as in RY Scuti ([Grundstrom et ahl 
12007ft . 

Dynamical interactions between a binary system and 
other stars can decrease the binary orbital period 
through angular m omentum tra nsfer, or through three- 
body interactions (Kozai 1962). Indeed the incidence 
of contact binary systems is hig her in stellar clusters 
than in the field (|Rucinskil 11998ft due to the more fre- 
quent three-body interactions. Recent models predict a 
timescale of 6-13 Gyr for an isolated low-mass binary 
system to undergo enough angular momentum loss to 



evolve from a s table orbital p eriod of ~l-2 days to be- 
low 0.22 days (Stepien 2006a. b). Estimates of angular 
momentum loss from single stars, applied to short pe- 
riod binary systems, show that magnetic braking takes 
on the order of a Hubble time to compress these peri- 
ods from 1-2 days to below the theoretical limit (e.g. 
iMaceroni fc Montalbanl 12004ft . With their discovery of 
sever al low mass syste ms below the predicted 0.22 day 
limit . iNefs et al.l (|2012[ ) demonstrated that initial orbital 
periods must have been ~1 day in order for magnetic 
braking schemes to explain the short period systems. 
INefs et ail ([2012D also show that the orbital evolution- 
ary timescales and initial conditions for these rare short 
period binaries are still largely unconstrained, and that 
many possible mechanisms may be at work in forming 
such systems. Once the stars reach a contact scenario, 
unstable mass transfer via the filled Roche lobes, and 
mass loss through the outer Lagrang e points, predicts a 
rapid evolution towards coalescence (iJiang et al.|[ 2011'). 
This type of over-contact evolution ([Rucinskil 11997ft is 
thought to be a formation path for blue straggler stars in 
stella r clusters ([Bradstreet fc Gu inan 1994Y lMateo et al] 
11990ft . and signatures include mass loss, changes in lumi- 
nosity, and orbital decay. 

Our measured period decay of approximately |P| ~ 8 
sec each year is larger than found for typical W UMa 
binary systems. The largest amplitude period decay 
measured for a contact binary is that of V1309 Sco, a 
merger system with an exponentially decreasing period 
with an amplitude of |P| ~300 sec year -1 ([Tylenda et al.l 
12011ft . This system exhibited an extreme nova-like out- 
burst, similar to V838 Mon, as it underwent a catas- 
trophic merging event in 2008. Assuming a linear period 
decay of |P| ~ 8 seconds per year, we would expect SDSS 
J001641-000925 to have a lifetime of -10 3 years. How- 
ever, if the period evolution follows an exponential de- 
cay as in the higher mass V1309 Sco system, the binary 
merger timescale may be as short as ~ 10 2 years. 

This amplitude of angular momentum loss required by 
the observed period decay can't be explained solely by 
magnetic braking, and the stars are too widely sepa- 
rated for_gravitational waves to expel significant energy 
(| C haul \l9 78ft . The angular momentum loss is likely due 
to the system being in an over-contact configuration and 
having filled its critical Roche surface, with mass be- 
ing expelled from the outer Lagrangian points. This 
would carry much more angular momentum from the sys- 
tem than magnetic breaking or winds alone, a nd would 
lead to a rapid coalescence into a single star (|Li et al.l 
2004). However, we cannot exclude the possibility of 
mass transfer from the primary to the secondary star 
with little to no loss of angular momentum to the system. 
This mass redistribution could ins t ead lea d to a binary 
evolution as described bv iStepienl (|2006aft . with the or- 
bital period increasing, and binary coalescence proceed- 
ing on a much more gradual nuclear or thermal timescale 
(|Li et al.l 12004ft . Our follow up study of the period evo- 
lution will help discern between these two possibilities. 

We emphasize the rarity of this system, as SDSS 
J001641-000925 is the first true over-contact M dwarf 
binary known. The critical question remaining for the 
system is whether it is evolving rapidly towards coales- 
cence due to the over-contact configuration. Additional 
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high precision follow-up photometry, as well as search- 
ing existing time domain surveys, is underway to con- 
clusively determine whether the period is decaying at an 
exponential rate, as was the case for V1309 Sco before 
its dramatic merger. 
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